Intermediate-mass fragments from fission and multifragmentation 

in the spallation of 136 Xe 
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A correlation between the production and kinematic properties of the fragments issued of fission 
and multifragmentation is established in the study of the reaction 136 Xe + hydrogen at 1 GeV 
per nucleon, measured in inverse kinematics at the FRagment Separator (GSI, Darmstadt). Such 
observables are analysed in a comprehensive study, selected as a function of the decay mode, and 
related to the isotopic properties of the fragments in the intermediate-mass region. Valuable in- 
formation can be deduced on the characteristics of the heaviest product in the reaction, which has 
been considered a fundamental observable for tagging the thermodynamic properties of finite nuclear 
systems. 



I. INTRODUCTION 

We present a comprehensive experimental study of 
the decay processes occurring at the multifragmentation 
threshold in the spallation of 136 Xe induced by 1 GeV 
protons, through inclusive kinematic and production ob- 
servables. In particular, we develop a selection method to 
identify and quantify the relative shares of the decay pro- 
cesses, fission and multifragmentation, which lead to the 
production of intermediate-mass fragments. Once estab- 
lished a correlation between production and kinematic 
quantities as a function of the isotopic observable and 
the decay mode, we analysed the characteristics of the 
intermediate-mass fragments issued of fission and mul- 
tifragmentation, with the further aim of determining a 
set of observables to benchmark and constrain elaborate 
reaction models. 

In rclativistic nuclear reactions, especially when in- 
duced by hadrons or light nuclei, kinematic and ther- 
mal observables probe the transition between different 
reaction mechanisms as a function of the bombarding 
energy Already some decades ago, several exper- 

iments were focused on spallation with the purpose of 
studying the reaction mechanism (for a non-exhaustive 
series of pioneering experiments see refs. 0, 0, H, 0, H, 
H, OH EH ) • In nucleon- induced reactions, at the incident 
energy of few hundred MeV per nucleon, an excited and 
fully equilibrated compound nucleus is formed, which de- 
cays by the emission of fragments in the form of particle 
evaporation and fission, according to the number of final 
states available for different mass splits [rH Il3j . At the 
incident energy of several GeV per nuclcon, the highly ex- 
cited composite nuclear system has also a large probabil- 
ity of decaying by multifragmentation; see ref . 11411 f or a 
large collection of reviews on this process, refs. [I5l.ll6l.ll7j 
for specific reviews treating multifragmentation induced 
by high-energy protons. In nucleon-induced reactions at 
incident energies which are intermediate between these 
extreme situations, for instance around 1 GeV per nu- 
cleon, the multifragmentation threshold could be at- 
tained: in this case, the exit channel of the system is 



determined by the competition between multifragmenta- 
tion and sequential decay. 

Such a complex phenomenology inspired the formula- 
tion of several alternative descriptions, attempted since 
the early eighties, in order to understand the transi- 
tion between the different reaction mechanisms. Ini- 
tially, a description of multifragmentation was proposed 
as a process of successive binary splits [18| within the 
transition-state model, or as an extension of the saddle- 
point interpretation to highly excited systems [l^ | . Later, 
several experimental data gave arguments for the large 
success of statistical models, within the freeze-out pic- 
ture [II HIL HI, [Hill, so that binary splits might even 
appear among the possible fragment configurations. Not 
only light-ion induced reactions comply largely with the 
statistical nature of the process (ISIS experiments [HI]), 
but also recent heavy-ion experiments at Fermi energy 
(INDRA [H, [M H3, HI US HI, MULTICS [H) gave 
special prominence to such a picture. As a function of 
the violence of the collision, the bimodal distribution of 
the largest-fragment size is interpreted as a signal of the 
finite-size counterpart of the liquid-gas phase transition 
in nuclear matte r |32[ [33|, which is associated to a fi- 
nite latent heat [3J]: the energy (and density) of the 
system and the size of the largest fragments are order 
parameters which rule the evolution of the system from 
a configuration where a large cluster survives at low ex- 
citation energy to a configuration where the system dis- 
assembles in small fragments of comparable size at high- 
excitation [H, HU . It should be pointed out that a dedi- 
cated study of the effect of the Coulomb field on the frag- 
mentation process further enriched this picture and re- 
vived the interest on the phenomenology of binary fission 
in highly excited systems. The phenomenon of Coulomb 
frustration, already well known in other systems, rang- 
ing from magnets on specific lattices to liquid crystals, 
from spin glasses to protein folding, can not minimise si- 
multaneously the long-range repulsive Coulomb interac- 
tion and the short-range attractive nuclear interaction, 
and results in one additional phase transition process, 
from the liquid phase to fission. This phenomenon of 
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Liquid-Fission transition is not associated to any latent 
heat. It was observed and detailed within the Liquid Gas 
Model [13, [H] as well as in phenomenological nuclear sta- 
tistical models [13, H^] • 

The experimental study detailed therein analyses the 
kinematic, production and isotopic observables of the two 
decay modes, fission and multifragmentation, tagged by 
the production of a specific nuclide in the intermediate- 
mass region. In the case of fission, such nuclide is the 
lightest binary partner, the other being close to the mass 
of the projectile. In the case of multifragmentation, such 
nuclide is one of the several fragments of similar size pro- 
duced by the division of the source and has only a fraction 
of the mass of the projectile. In both cases, indirectly or 
more directly, the experimental signature we obtain by 
the kinematics of the intermediate-mass fragments gives 
a good information on the size of the heaviest product in 
the reaction. Especially in the case of multifragmenta- 
tion, the characteristics of the heaviest fragment have a 
fundamental role in carrying robust signals of the finite- 
size counterpart of the nuclear liquid-gas phase transi- 
tion. From such premises, the main aim of this experi- 
mental work is threefold: 

1. Analysis. We focus on the selection problematics 
for high-resolution inclusive experiments, and we 
propose a new analysis procedure for correlating 
production and kinematic observables. 

2. Experimental results. We provide several inclusive 
intermediate-mass fragment observables, MS i\ func- 
tion of the decay process: we quantify both isotopic 
cross sections and invariant cross sections, concern- 
ing the production and the kinematics, respectively. 

3. Discussion. Concerning the exit channel: we anal- 
yse the boost imparted to the partners of binary 
splits, we extract the distribution of fragments 
along the mass number, and we analyse the strong 
staggering which characterises the nuclide produc- 
tion. Concerning the connection between the en- 
trance channel and the exit channel: from the kine- 
matic observable, we propose a new indicator for 
the impact parameter and we illustrate a possible 
interpretation for the evolution of the kinematic ob- 
servable MS M function of the size of the fragments 
for the different decay modes. 



II. 



EXPERIMENTAL OVERVIEW 



The contribution of two processes to the decay, fis- 
sion and multifragmentation, makes experimental con- 
ditions rather challenging. Inclusive experimental ap- 
proaches, based on the use of magnetic spectrometers 
are specially adapted to relativistic nuclear reactions in- 
duced by hadrons: the fraction of the bombarding en- 
ergy which remains in the system appears almost com- 
pletely as thermal energy, and the available phase space 



is fully explored. These clear experimental conditions, 
reflected in the isotropic kinematics, make the inclusive 
approach suitable for connecting the production observ- 
ables to the kinematic observables: in particular, the 
production observables consist in the full isotopic iden- 
tification, extended to the complete range of fragment 
sizes; the kinematic observables are the distribution of 
longitudinal velocities, measured for single nuclides, from 
the magnetic rigidity, with very high precision and with- 
out any threshold. Nevertheless, we show in the follow- 
ing that some difficulties still remain when rather heavy 
fragments are studied, because the analysis procedure we 
adopted, based on the disentangling of the full velocity 
distribution in more components, becomes limited. 



A. Spallation at the multifragmentation threshold 

A first inclusive technique of disentangling reaction 
processes on the basis of longitudinal-velocity spectra 
was introduced for distinguishing the average character- 
istics of the fission fragments from those of the evapora- 
tion residues [i^ . [5fl |. in inverse-kinematic experiments 
at the FRagment Separator [42|, [H[ at GSI in Darm- 
stadt. Concerning spallation at the multifragmentation 
threshold, still with the same facility and experimen- 
tal approach, the kinematic and production properties 
of intermediate-mass fragments were deduced in an in- 
clusive experiment for the reaction 56 Fe+ hydrogen at 1 
GeV per nucleon [i(| With the purpose of recon- 
structing the invariant-velocity spectra of intermediate- 
mass fragments from inclusive kinematic observables, a 
simple procedure of analysis of the longitudinal-velocity 
spectra was introduced: it was based on the assump- 
tion that all exit channels were associated to the same 
value of the linear momentum transfer for one individ- 
ual nuclide formed as a final product. The experimen- 
tal results revealed mixed signals of multifragmentation 
and fission [40| . which could however not be quantified 
in their relative share and identified separately; a pos- 
sible disentangling of the two exit channels in terms of 
kinematic and thermal properties could be proposed only 
through the comparison with a statistical model. The 
same system was also the object of an exclusive exper- 
iment at ALADIN 



44[ at GSI. Profiting from particle 



correlations, still not reinforced by kinematic observables, 
it confirmed the mostly binary character of the produc- 
tion of intermediate-mass fragments and the presence of 
a few decays with larger fragment multiplicity. Also in 
this case, the analysis of the reaction mechanism was dis- 
cussed on the basis of statistical models. 

In a more recent measurement at the FRagment Sep- 
arator, the same inclusive experimental approach previ- 
ously used for iron, was followed to measure the system 
136 Xe + hydrogen at 1 GeV per nucleon (H. 136 Xe 
is the stable nucleus with the largest neutron excess 
N — Z among the nuclei below the Businaro-Gallonc 
point [46|, , where the saddle point becomes unstable 
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FIG. 1: (Color online) Intra- nuclear-cascade calculation of 
the distribution of excitation energy of the hot fragments in 
the system 136 Xe(! a GcV) + P\ above around 3 MeV per nu- 
cleon the onset of multifragmentation is expected. Elsewhere, 
fission-evaporation is the dominant decay process. 

to asymmetric splits. Binary decay should result in the 
production of neutron-rich intermediate-mass fragments. 
A large feeding of evaporation-fission decays should re- 
sult from the range of excitation energies explored by 
this system. To corroborate our expectations, fig. [1] 
presents a calculation within the intranuclear exciton- 
cascade model [48| of the excitation distribution, which 
extends to large values: a remarkable fraction of the pro- 
duction extends to excitation energies exceeding 3 MeV 
per nuclcon. According to the large amount of theoreti- 
cal and experimental results, the transition towards mul- 
tifragmentation is expected above this value and, again, 
similarly to fission, it should become manifest in the for- 
mation of neutron-rich intermediate-mass fragments. For 
this reason, measurements of production cross sections 
alone (and in particular the observables related to the 
neutron enrichment), uncorrelated to the measurement 
of kinematic observables are neither adapted for probing 
the initial excitation process nor for selecting the exit 
channel in such kind of systems. Additional informa- 
tion on the reaction mechanism is necessary and it can 
be deduced either from complete correlations in exclu- 
sive experiments or, as we illustrate in the present work, 
from high-resolution velocity spectra of fragments in an 
inclusive experiment. 



B. Inclusive event selection 

We detail in this paper the new method of inclusive 
event selection which we developed for the analysis of 
the system 136 Xe+p. In comparison with the analysis 
method used for the system 56 Fe+p, the new method has 
the advantage of being a quantitative technique of event 
selection, allowing to associate the probability for a given 
decay pattern to a group of events. In particular, the 
inclusive event selection method is used to identify kine- 
matically the two processes of fission and multifragmen- 
tation in the formation of intermediate-mass fragments, 
and to quantify them in terms of nuclide cross sections 
and kinematic properties, profiting from the two main 



specificities of the experimental method we employed. 
The first specificity is the possibility to measure small 
velocities in the projectile frame with high resolution, so 
as to access the physics of the small-energy portions of 
the kinetic-energy spectra of fragments; this portion is 
found in this work to give constraints for the identifica- 
tion of the exit channel. The second specificity is the 
possibility to connect such kinematic observable to the 
production of a given nuclide. 



C. The FRagment Separator 

The fragments were produced in inverse kinematics by 
directing a primary beam of 136 Xe at 1 A GeV on a target 
of liquid hydrogen contained in a cryostat with thin tita- 
nium windows. The projectile residues were analysed in- 
flight, with the Fragment Separator [H,[43|]. Four dipole 
magnets compose the magnetic spectrometer, a disper- 
sive focal plane is localised behind the first two dipole 
magnets, and the full spectrometer was set achromatic. 
The identification of a fragment was deduced from the 
measurement of its time of flight, from recording the po- 
sition where its trajectory intersects the dispersive focal 
plane, and from its nuclear charge measured in ionisa- 
tion chambers, placed in the proximity of the terminal 
focal plane. Recoil velocities and kinetic energies were 
deduced from the magnetic rigidities, which were mea- 
sured by the spectrometer with a relative uncertainty of 
5 • 10~ 4 (FWHM) for individual reaction products. 

A series of corrections and calibrations were required 
for deducing the distribution of longitudinal recoil veloc- 
ities for each nuclide. The raw spectrum was slightly 
corrected by taking into account the energy loss in the 
target. The events related to reactions in any layer of 
matter which does not correspond to hydrogen were sup- 
pressed; these events come for the greatest part from the 
material (mostly titanium) constituting the cryostat, and 
their contribution was measured with a specific experi- 
ment. The number of counts was normalised to the beam 
dose per target thickness, which was deduced by measur- 
ing the current of the primary beam. The data treat- 
ment required at this point to correct for the decrease of 
the angular transmission for ions with magnetic rigidities 
which deviate from the nominal values. This effect is seen 
in the velocity spectra, when they are constructed from 
composing different magnetic settings of the spectrome- 
ter: they present structures which appear periodically for 
each portion of the spectrum corresponding to one mag- 
netic setting (40| . Since the angular transmission can be 
formulated simply as a function of the position of the 
trajectory intersectin g th e dispersive focal plane and the 
terminal focal plane [60(, these structures are removed 
by correcting the data in a way that they correspond to 
the transmission of ions with nominal magnetic rigidi- 
ties. Finally, additional corrections took into account 
the modification of the yields due to secondary reactions 
occurring in the hydrogen target and in other layers of 
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FIG. 2: (Color online) Systematics 
of measured velocity distributions as 
a function of the nuclide for elements 
ranging from Li to Si, indicated with 
black dots; error bars show the sta- 
tistical uncertainties. The measured 
spectra are disentangled by a fitting 
procedure (see section ITTi]) in the con- 
tribution of a coulomb-hole mode and 
a multifragmentation mode (see sec- 
tion II V|) . The shapes of these two 
components, which are concave and 
convex, respectively, are deduced so 
that their sum is a fit of the measured 
shape. 
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matter intercepting the beam of fragments. These correc- 
tions are discussed in details in ref. [45[ and they are the 
standard procedure of analysis for all nuclidc-production 
measurements at the FRagmcnt Separator. 

At this stage, a distribution of measured events as a 
function of the longitudinal recoil velocity in the pro- 
jectile frame is associated to each reaction product, 
identified in mass number and nuclear charge. The in- 
tegral I of each distribution is equal to the fraction of 
production cross section which is selected by the angular 
acceptance of the spectrometer. Fig.[3]shows a systemat- 



ics of measured longitudinal- velocity spectra dl(v^)/dv^ 
normalised to the same integral for all measured isotopes 
of the elements ranging from lithium to silicon. In or- 
der to obtain invariant-cross-section distributions inde- 
pendently of the experimental conditions, the geometry 
of the angular acceptance of the spectrometer needs to 
be accurately accounted for in the data analysis. 

In the following section, we describe the method for re- 
constructing the invariant-cross-section distribution from 
the measured longitudinal-velocity spectra. Such tech- 
nique is a further development of the method used in the 
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analysis of the spallation of iron at 1 GeV per nucleoli 
measured at the FRagment Separator [ipj . 



III. DECONVOLUTION OF THE INCLUSIVE 
VELOCITY SPECTRA 

When the longitudinal-velocity distribution recorded 
within the limited angular acceptance of the spectrome- 
ter is measured with very high resolution, valuable infor- 
mation on the reaction pattern can be recovered by its 
deconvolution. This procedure is based on the picture 
of the reaction process proceeding by a fast stage, char- 
acterised by nucleon-nucleon collisions, which ends up in 
an excited system, and a slower deexcitation phase [5l|. 
The interaction with the target proton induces a mean 
shift and a broadening of the velocity distribution of the 
projectile. The shift is expected to be a function of the 
impact parameter [l[ . The broadening is well represented 
by the Fermi momenta of the removed nucleons [52| and, 
within this picture, is isotropic [53|. In the course of the 
deexcitation, the system is considered thermalised. This 
justifies simplified assumptions on the angular distribu- 
tion of the emitted particles. After thcrmalisation, the 
angular distributions of the particles emitted in sequen- 
tial decay that reaches from evaporation of light particles 
up to symmetric fission, or produced by multifragmcnta- 
tion, are forward-backward symmetric in the system of 
the emitting source. The angular anisotropy in the emis- 
sion of light charged particles can be estimated within 
the Hauser-Feshbach picture [54|, the anisotropy of fis- 
sion processes can be estimated in the framework of the 
statistical population of the K-quantum numbers [55| and 
the anisotropic character of the multifragmcntation pro- 
cess was already successfully described in early transport 
models l57l [58j. However, since the angular momen- 
tum introduced in the collision phase is rather low and 
the excitation energy rather large, for the deconvolution 
procedure, we assume the angular distribution to be ap- 
proximately isotropic for any decay process, in the frame 
of the decaying source. In agreement with this expecta- 
tion, we add that the reaction analysed in this work is in 
between a large range of incident energies for which signa- 
tures of isotropy were measured in several experiments: 
from lowly excited fissioning systems [l3| to spallation 
with high-energy light projectiles and relativistic heavy- 
ion collisions [lg, llZl, [59|)- Still, we anticipate that we 
will find slight deviations from the isotropy in the data, 
which can not be described within our assumptions. The 
deconvolution procedure is sketched for the nuclide 20 F 
in fig. [3] and is essentially based on three points: 

1. Especially when dealing with relativistic protons, 
the kinematics of one individual nuclide can be 
decomposed in several isotropic emission patterns, 
each one centred around one apparent source u. 

2. Each apparent source u is localised at a longitudi- 
nal velocity vP in the projectile frame (and u L in 



the laboratory frame), which reflects the momen- 
tum transfer involved in the collision and is aligned 
along the beam axis: this assumption is equivalent 
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FIG. 3: (Color online) Procedure of reconstruction of the 
invariant cross sections applied to the nuclide 20 F. The de- 
composition in the concave component of the Coulomb-hole 
mode and the convex component of the multifragmentation 
mode indicated by different line types as in fig. [2] (a) Mea- 
sured velocity distributions (black dots) and fitted spectra as 
in fig. [2] The decomposition of the multifragmentation mode 
in three Gaussian components is shown, (b) Planar cut along 
the beam axis of the reconstructed full distribution d<r/di; p 
in the beam frame as a cluster plot. The two lines indicate 
the boundaries of the angular acceptance, inside of which the 
fragments could be measured. The ridge of the Coulomb- 
hole component is indicated by the circle of radius i> pca k- (c) 
Projection on the beam axis of the reconstructed full distri- 
bution da/dv p . Arrows indicate the positions of the source of 
the Coulomb-hole component u p , the average source position 
of the multifragmentation component (u^ ) and the average 
longitudinal velocity (vp. (d) Invariant-cross-section distri- 
bution, (e) Reconstructed cross-section distribution for the 
Coulomb-hole component as a function of the boost velocity 
in the source frame Dboost- The plot allows to deduce the 
mean boost (uboost), which differs from w pea k and from « rm8 . 
(f) Reconstructed cross-section distribution for the Coulomb- 
hole component as a function of the total kinetic energy in 
the source frame Eboost- v ruls is deduced from the mean value 

(-Eboost)- 
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to imposing that the overall kinematics is symmet- 
ric with respect to the beam axis. 

3. The shape of the angular acceptance is known pre- 
cisely. 

The first two points are assumptions. The third point is 
an ingredient guaranteed by ion-optical calculations [6(| ■ 
As described in section HHthese ion-optical calculations 
were also included in the procedure of raw-data treat- 
ment which had as an outcome the measured spectra of 
fig. [U they are validated by requiring that the different 
magnetic settings of the spectrometer do not produce pe- 
riodic structures on these spectra (4fjj ; such a requirement 
is a tool to check and even refine the exact properties of 
the angular acceptance. 

We already defined the quantity d/(v p )/dv p as the 
measured longitudinal-velocity spectrum in the projec- 
tile frame (see fig. [2|; it is a sum of the contributions 
of sources complying to the hypotheses (f) and (2), and 
is affected by the angular-acceptance shape of the spec- 
trometer. Additionally, for a given number of sources of 
velocity u P in the projectile frame, since the velocities in 
the source frame are not relativistic, we define 



da 



d(v. 



u 



') 



(1) 



as the invariant cross section attributed to the source i. 
This is a production quantity which does not depend on 
the angular acceptance and, according to the assumption 
(1), it depends only on the absolute value of the relative 
velocity u p — m p in the frame of the source u P . 

These cross sections are obtained by the deconvolu- 
tion of the measured spectrum d/(i> p )/di> p taking into 
account the ion-optics required by the point (3). This 
procedure is detailed in a general form in appendix [A] 
By composing these contributions, we construct the 
invariant- velocity cross section accounting for all sources, 
independently of the angular acceptance, which is de- 
noted by 



oi(t; p ) 



c 2 da 



As explained in appendix [AJ we obtain the relation 



(2) 



(3) 



The spectrum reduces to a fully isotropic distribution 
if all emitting sources coincide for one individual nuclide. 

This prescription was introduced in the simplified form 
of isotropic emission from one unique source in the analy- 
sis of the light residues in the system 56 Fe(! a GcV) +P EOl • 
For such system, no evidence was found for postulating 
a spreading of emitting sources. On the contrary, the 
description of multiple-source emission was necessary in 
the analysis of intermediate-mass fragments produced in 
the spallation of xenon. 



IV. RESULTS 

As shown in fig. [2l the evolution of the measured 
longitudinal-velocity spectra over a systcmatics of sev- 
eral light nuclides ranging from lithium to silicon invites 
to disentangle two main components: one has a shape 
with a concave centre boarded by two maxima: it cor- 
responds to a Coulomb hole; the other component has a 
convex shape with a tail extending in the backward di- 
rection: as discussed in the following, we attribute to this 
mode the features of a multifragmentation process. The 
decomposition of the measured longitudinal- velocity dis- 
tribution in the Coulomb-hole mode (symbol ^) and the 
multifragmentation mode (symbol ^) is shown in fig. [2] 
for elements ranging from lithium to silicon. The two 
kinematic modes are disentangleg by a fitting procedure 
which respects the constraints imposed by the detailed 
knowledge of the ion-optics. In the plot (a) of fig. [3] 
this decomposition is shown for 20 F. The convex shape of 
the multifragmentation mode is fitted as a skewed Gaus- 
sian distribution and the concave shape of the Coulomb- 
hole mode is fitted as the superposition of two Gaussian 
shapes with integrals constrained to have a ratio com- 
patible with the effect of the angular acceptance on the 
forward and backward emission kinematics. The forward 
side of the spectra of the heavier isotopes of beryllium 
were not fully measured or partially affected because the 
magnetic rigidity of the primary beam was circumvented 
in the experiment; for these nuclides, the constraint on 
the concave shape of the Coulomb-hole mode allowed to 
profit from the precise measurement of the backward side, 
and the parameters for the convex shape of the multifrag- 
mentation mode were deduced from isobaric nuclides. 

The double-humped shape of the Coulomb-hole com- 
ponent, with integral i^(i> p ), was interpreted as a shell in 
velocity space centred around one single isotropic emit- 
ting source of recoil velocity u p in the beam frame. The 
plot (b) of fig. [3] shows a planar cut along the beam di- 
rection of the reconstructed velocity distribution da/dv p ; 
the radius v pea k indicated in the figure is the radius of the 
ridge of the ring marking the Coulomb- hole component. 
The same figure shows that the distribution da/dv p is 
cut by the conical boundary of the acceptance so as to 
produce the two humps in the longitudinal- velocity spec- 
tra. The larger is the radius of the emission shell, the 
larger is the integral of the forward hump with respect 
to the backward hump. As the two peaks are spaced of a 
distance which largely exceeds their widths, their shapes 
are close to Gaussian distributions. 



Invariant cross sections 



The asymmetric shape of the multifragmentation com- 
ponent, of integral J^(w p ), could not be interpreted as 
the emission from one single isotropic source: we decom- 
posed it in more Gaussian components dli^,(v^)/dv^ of 
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FIG. 4: (Color online) Systematics of 
invariant-cross-section distributions as 
a function of the nuclide for elements 
ranging from Li to Si decomposed in 
the Coulomb-hole and multifragmen- 
tation components. 
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different integrals, and different mean values uP ; . We 

assumed that the emission kinematics is independent of 
the recoil velocity of the corresponding source u p : this 
assumption results in assigning to all the Gaussian com- 
ponents the same width, which is determined by the for- 
ward side of the convex shape of the multifragmentation 
mode (the effect of the angular acceptance on a Gaussian- 
like emission kinematics is discussed in ref. [(3(j). In the 
plot (a) of fig. [31 a decomposition of the multifragmen- 
tation mode in three Gaussian components is shown; for 
this spectrum, a number of three Gaussians was chosen so 
that their spacing is smaller than their width, sufficiently 



to avoid the appearing of structures. 

By the fitting procedure, the measured spectrum is dis- 
entangled in the Coulomb-hole and multifragmentation 
components, so that 

dv* duP ^ Av\ ' {) 

I II i 

All parameters involved in the fit have uncertainties de- 
duced from the corresponding covariance matrix. The 
plot (b) of fig. [3] confirms that the backward-forward 
asymmetry which characterises the Coulomb-hole com- 
ponent in the measured spectra is produced by the angu- 
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lar acceptance of the spectrometer and disappears in the 
reconstructed distribution. On the contrary, the asym- 
metry of the multifragmcntation component is not an 
effect of the acceptance and is kept in the reconstructed 
shape: it depends therefore on the reaction process. 

The projection of the reconstructed distribution 
da/dv p on the beam axis is illustrated in the plot (c) 
of fig. 02 it corresponds to the longitudinal- velocity dis- 
tribution we would measure if the spectrometer had a 
full acceptance. In this representation the integrals of 
the different kinematic components give the correspond- 
ing production cross sections and their mean values give 
the average recoil of the corresponding sources. 

The variation of the reconstructed quantity da/dv p 
along the beam axis coincides with the distribution of the 
invariant cross section cri (^[j*) illustrated in the plot (d) of 
fig. [3l The Coulomb-hole component <ti ^ and the multi- 
fragmentation components a\ . ^ of the distribution of in- 
variant cross section are deduced from the corresponding 
components of the measured longitudinal-velocity spec- 
tra. The procedure of dcconvolution requires the knowl- 
edge of the mean velocities of the sources x£ and u p ^. 
These quantities are however not measured directly, but 
deduced through an optimisation procedure. The com- 
plete spectrum of invariant cross section is reconstructed 
by summing its components, shifted of the corresponding 
positions of the sources: 



CI. 



o + 



a, (5) 



This is the application of Eq. (02) when v p is constrained 
to vary along the beam axis. The procedure of velocity 
reconstruction was applied for all intermediate-mass frag- 
ments ranging from lithium to silicon. In fig.|Uthe result- 
ing systematics of invariant-cross-section distributions is 
shown, deduced from the deconvolution of the measured 
longitudinal-velocity spectra illustrated in fig. 03 The 
main parameters characterising the reconstructed spec- 
tra are listed in table IIIII in the appendix [B] 

One more representation of the Coulomb-hole compo- 
nent is obtained by integrating the corresponding recon- 
structed distribution da/dv p over the polar angle. This 
is shown in the plot (e) of fig. El it probes the evolu- 
tion of the cross section as a function of the Coulomb 
boost Wboost- The plot (f) of fig. [3] is an equivalent rep- 
resentation in terms of the total kinetic energy Eb oost . 
The mean value (-©boost) is indicated in the plot (f) and 
corresponds to the velocity v rms indicated in the plot (e). 
v Ims = y/ (wboost) 2 + V accounts for the variance V of the 
velocity distribution and is slightly larger than the mean 
velocity (wboost), also indicated in the plot (e). In section 
IV Al the quantity v rms is compared to fission velocities 
calculated under different prescriptions for the Coulomb 
boost. 



B. Production cross sections 

The production cross sections were deduced from in- 
tegrating the reconstructed distributions da/dv p . The 
disentangling of the velocity spectra into the two kine- 
matic modes allowed to separate the contributions of each 
mode. The results are presented for the isotopic chains 
of elements ranging from lithium to silicon in fig. [5] and 
listed in table |TJ The uncertainties account for statisti- 
cal and systematical contributions (see ref. [45| for de- 
tails). The production of the lightest elements is mostly 
related to the Coulomb-hole mode. For heavier elements 
the multifragmcntation mode gains gradually a larger 
fraction of cross section and it equals the Coulomb-hole 
mode around silicon, where the method of deconvolution 
reaches its limit of applicability. This evolution can be 
appreciated in the projections of the cross sections along 
the element and neutron number, in fig. [5] The average 
neutron-to-proton ratio (N/Z) of the intermediate-mass 
fragments produced by the two modes are indistinguish- 
able, as shown in fig. [Jj both the two modes feed preva- 
lently the neutron-rich side of the nuclide chart with very 
similar shapes of the cross-section distributions. 
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FIG. 5: (Color online) Production cross sections for the iso- 
topes of elements ranging from Li to Si. The contributions 
from the two kinematic modes are shown. 
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TABLE I: Nuclide cross sections a measured in this work for the production of elements ranging from Li to Si. 
section restricted to the Coulomb-hole or multifragmentation mode is indicated as or , respectively. 



The cross 



A a [mb] [mb] a^. [mb] A a [mb] 



[mbl 



[mbl 



A a [mb] 



[mb] 



[mb] 



9 
10 
11 
12 

10 
11 
12 
13 

11 
12 
13 
14 
15 

13 
14 



Li 

10.12(59) 

19.0 ± 1.1 

3.90(30) 
Be 

4.20(25) 
4.61(28) 
0.290(23) 
0.142(10) 
B 

1.131(73) 
4.88(29) 
1.33(10) 
0.406(27) 
C 

0.212(28) 
1.528(92) 
2.15(13) 
1.65(12) 
0.176(27) 
N 

0.0241(28) 
0.446(33) 



9.36(56) 0.757(60) 

18.0 ± 1.1 1.066(85) 

3.64(29) 0.259(33) 

3.82(23) 0.378(32) 

4.22(26) 0.392(32) 

0.220(19) 0.070(10) 

0.1267(94) 0.0151(20) 

0.999(66) 0.132(18) 

4.25(26) 0.638(46) 

1.055(89) 0.277(32) 

0.355(25) 0.0512(55) 

0.166(26) 0.0459(75) 

1.350(84) 0.178(13) 

1.90(12) 0.247(27) 

1.197(89) 0.449(52) 

0.135(25) 0.0412(80) 

0.0197(27) 0.00436(59) 

0.374(30) 0.0723(62) 



15 
16 
17 

15 
16 
17 
18 
19 

17 

18 
19 
20 
21 

19 
20 
21 
22 
23 

22 



1.85(13) 
0.499(38) 
0.276(50) 
O 

0.039(13) 
0.603(40) 
0.534(38) 
0.656(50) 
0.256(33) 
F 

0.00925(79) 
0.095(10) 
0.388(28) 
0.494(44) 
0.404(70) 
Ne 

0.0068(24) 
0.126(13) 
0.336(28) 
0.538(44) 
0.247(45) 
Na 

0.0677(80) 



1.42(11) 

0.416(35) 

0.228(48) 

0.032(13) 
0.485(35) 
0.350(28) 
0.525(45) 
0.168(30) 

0.00376(23) 

0.0678(91) 

0.281(23) 

0.358(39) 

0.285(65) 

0.0049(24) 

0.078(12) 

0.223(21) 

0.392(36) 

0.189(43) 



0.422(36) 

0.0823(87) 

0.048(11) 

0.0075(18) 

0.118(10) 

0.184(19) 

0.131(12) 

0.088(11) 

0.00549(69) 

0.0266(28) 

0.1066(93) 

0.136(13) 

0.119(22) 

0.00186(47) 

0.0479(46) 

0.113(14) 

0.146(19) 

0.058(11) 



0.0480(71) 0.0197(28) 



23 0.321(22) 

24 0.305(21) 
0.282(35) 
Mg 

0.00559(46) 
0.1196(91) 
0.271(27) 
0.394(44) 
27 0.242(72) 
Al 

0.0015(38) 
0.0453(56) 
0.271(29) 
0.294(34) 
0.268(31) 
0.092(13) 
Si 

0.080(11) 
0.199(27) 
0.427(53) 
0.199(22) 
0.113(16) 



25 

23 
24 
25 
26 



25 
26 
27 
28 
29 
30 

28 
29 
30 
31 
32 



0.205(13) 0.116(14) 

0.191(12) 0.114(14) 

0.208(32) 0.074(10) 

0.00384(29) 0.00175(29) 

0.0723(46) 0.0473(65) 



0.1400(87) 

0.232(20) 

0.132(62) 

0.0012(38) 

0.0255(22) 

0.159(14) 

0.164(14) 

0.153(14) 

0.0463(38) 



0.131(23) 
0.162(36) 
0.110(36) 

0.00032(63) 

0.0197(49) 

0.112(24) 

0.130(29) 

0.115(26) 

0.046(12) 



0.0388(33) 0.042(10) 

0.0977(81) 0.101(25) 

0.233(21) 0.194(46) 

0.123(10) 0.076(18) 

0.0555(47) 0.057(14) 



C. Limitations of the method 

It was not possible to extend the deconvolution for 
heavier elements than silicon because the uncertainties 
of the fitting parameters for the identification of the 
two kinematic modes in the measured velocity spectra 
determined a limit for the analysis. Above silicon, a 
simplified method was employed, based on the assump- 
tion of one unique isotropic source associated to the 
kinematics of one fragment [45[ ; the same simplification 
was also employed in the analysis of the kinematics of 
the intermediate-mass fragments formed in the system 
56 Fe (1AGcV) +pM- 

Concerning the multifragmentation mode, we already 
pointed out the presence of a tail extending to the back- 
ward direction which characterises the lightest fragments. 
Such tail should result from the early stages of the col- 
lision and reflects fluctuations in the momentum trans- 
fer: within the intranuclear cascade model, these fluc- 
tuations are connected to a smaller or larger interaction 
of the shower of ejectile nucleons with the system, and 
result in a distribution of emitting sources. This distribu- 
tion would be attributed to the longitudinal momentum 
transfer due to friction, while the shape of the isotropic 
distribution would be governed by the Fermi momenta 
of the emitted nucleons [12] and the recoil in the de- 
cay process (fission or multifragmentation). It must be 
admitted, however, that the relatively small angular ac- 
ceptance does not allow to get direct experimental infor- 



mation on the shape of the velocity distribution in trans- 
verse direction. Guided by the present knowledge from 
other experiments and in agreement with the theoretical 
expectations for the contributions from Fermi momenta 
and decay recoil, we still keep the underlying assump- 
tion of isotropic emission patterns as a technical option 
to describe the tail to backward velocities. However, the 
validity of this assumption may only be tested by a full- 
acceptance experiment. 



V. DISCUSSION 

In the analysis of the velocity spectra of the 
intermediate-mass fragments formed in the system 
56 Fe(! a GcV) +P E3| similar results were found: the effect 
of the Coulomb field on the kinematics of fragments re- 
veals the interplay of a Coulomb-hole and multifragmen- 
tation mode, as shown in the left column of fig. [HI How- 
ever, in such system the Coulomb-hole mode exhibits two 
humps with large widths, with the consequence that the 
two modes could not be disentangled in single spectra and 
they were identified only for the extreme case when the 
contribution of one component is dominant; their coex- 
istence in the production of one individual nuclide could 
only be postulated on the basis of a model calculation, 
as shown in the right column of fig. [H The calculation, 
where the distribution of hot fragments is evaluated with 
the intranuclear exciton-cascadc model [48[, and their 
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FIG. 6: (Color online) Projections of the intermediate-mass- 
fragment cross sections along the element and neutron num- 
ber. The contribution from the two kinematic modes are in- 
dicated. 



1.5 
1.4 
1.3 
1.2 
1.1 
1 



136 Xe 



o Multifrag. 

Coulomb 
hole 

x Sum 








12 16 20 24 28 
Z 



FIG. 7: (Color online) Evolution of the average neutron-to- 
proton ratio (N/Z) of the intermediate-mass fragments as a 
function of the element number. The contributions from the 
two kinematic modes are indicated. 



breakup is simulated by the Statistical Multifragmen- 
tation Model [2(| HH HI!, associates the Coulomb-hole 
component in the velocity spectrum to events with mul- 
tiplicity of larger fragments than alpha M — 2 and the 
multifragmcntation component to events with multiplic- 
ity M > 2. 

On the contrary for the heavier system 136 Xe(i a Gov) + 
p, the identification of the two distinct channels in the 
formation of intermediate-mass fragments is an experi- 
mental result, even when they coexist in the formation 



of the same nuclide. To our knowledge, there have not 
been previous attempts to extract this observable from 
inclusive experiments and for similar systems. 



A. Exit channel 

1. Coulomb hole 

Double-humped spectra were measured at the FRag- 
ment Separator in several experiments: they were usu- 
ally taken as the signal of binary fission characterised by 
mostly symmetric splits (without being exhaustive, we 
cite refs. [H, and were also observed in presence of 
asymmetric fission splits, when light fragments were mea- 
sured [6l| . With the purpose of comparing to fission data 
it is therefore interesting to quantify the Coulomb repul- 
sion reflected by the forward and backward peaks in the 
present data. In fig. [SI the boost v rms is deduced as in the 
plot (e) of fig. [3] and its evolution is studied as a function 
of the neutron number for elements ranging from lithium 
to silicon. The same experimental points are collected 
in table HT1 As described in section HV A[ v ims is the ve- 
locity deduced from the mean value of the total kinetic 
energy distribution (-Eboost), which represents the aver- 
age mechanism leading to the observed nucleus. Because 
of the wide velocity spectra, this quantity results slightly 
larger than the mean value of the velocity distribution 
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FIG. 8: (Color online) Left column. Measured velocity distri- 
butions and invariant-cross-section distributions for four nu- 
clides produced in 56 Fe(i a GeV) +P- Right column. Invariant- 
cross-section distribution for 6 Li, calculated with the intranu- 
clear exciton-cascade model l4 p|. c oupled with the Statistical 
Multifragmentation Model [201. 12U [23| ; two components, as- 
sociated to the multiplicity M = 2 and M > 2 (for larger 
fragments than alpha) are indicated. Data and calculation 
from [H. 
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(vboost), which is also plotted in fig. [9] for comparison. 
In the figure, the reconstructed boost v lms is compared 
to the fission velocity calculated under different prescrip- 
tions, all applied to the same extreme case: the binary 
split of the largest possible mother nucleus 136 Xe, which 
also determines the largest boost. The calculations do 
not change sensibly for the split of lighter isotopes of 
xenon, if fission occurred after the emission of some neu- 
trons. Three prescriptions are tested. First, we com- 
pare to the empirical systematics of total kinetic energy 
for light fissioning nuclei condensed in the parametrisa- 
tion of Tavares and Terranova [62| . and further modi- 
fied in Ref. to describe asymmetric splits. Then, 
we compare to the scission-point model of Wilkins et al. , 
which was already used to describe similar data [HH . The 
model deduces the total kinetic energy from an empirical 
liquid-drop description, based on the parametrisation of 
Refs. [63|,[64||. The underlying picture of the model, which 
is the fission of a deformed system through the formation 
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FIG. 9: (Color online) Evolution of the reconstructed boost 
Dims of the Coulomb-hole mode (and of the mean boost 
(utmost), plotted for comparison, without error bars) as a 
function of the neutron number, for elements ranging from 
lithium to silicon. Three prescriptions for calculating the 
Coulomb boost in the binary split of 136 Xe are compared 
to the measurement: the total-kinetic-energy systematics 
of Tavares and Terranova [g^), the scission-point model of 
Wilkins et al. [6^, [64| and the nucleus-nucleus fusion model 

of Bass [Hi lei. 



TABLE II: Reconstructed boost ii rms , deduced for the 
Coulomb-hole mode. The uncertainty, indicated in paren- 
thesis, includes both statistical and systematical errors. 
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of a neck, is appropriate to describe the symmetric split 
of heavy fissioning systems, but is not sufficiently justi- 
fied when the formation of intermediate-mass fragments 
is described. This second calculation may be regarded 
as a lower limit. As a third option, we deduce the boost 
from the nuclear potential which the light fragment and 
its heavy partner would explore in a fusion reaction, ac- 
cording to the empirical model of Bass [65|, : within 
this more appropriate description, the Coulomb repul- 
sion acts on non-deformed fission fragments and results 
in the largest boost, compared to the first two prescrip- 
tions. For the extreme case of a fissioning 136 Xe, while 
isotopes from nitrogen to neon are well described both by 
the fission systematics and the Bass potential, the heav- 
ier isotopes from sodium to silicon could be compatible 
with a lighter fissioning system, according to the same 
two prescriptions, or still with the split of 136 Xe within 
the scission-point model of Wilkins et al. The large in- 
crease of the boost v rms for decreasing element numbers 
in the region from lithium to carbon, is not even de- 
scribed by the Bass potential. This phenomenology is 
rather general: the same conclusion was presented for 
the analysis of the system 56 Fe( lj 4 GcV ) +p and a 
similar tendency was found for the fission fragments in 
the system 238 U(i,4GoV) +P [HTJ] . Even the simultane- 
ous emission of a third light fragment or particle would 
not deform the shape of the Coulomb-hole component 
sensibly. The large value of the Coulomb boost v rms 
may manifest the presence of some additional contribu- 
tion to the total kinetic energy. To give a more general 
description, we relate the Coulomb-hole mode to events 
with low fragment multiplicity and large mass asymmetry 
(Ai — A2)(Ai +A2), where A\ is the mass number of the 
heaviest fragment and Ai indicates the second heaviest 
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fragment, the one the velocity spectrum belongs to. The 
large width of the backward and forward peaks {^ poak }: 
plotted in fig. [TU] for A < 30, is compatible with a large 
number of decay paths all ending in the production of the 
observed nuclide; these decay paths account for several 
different breakup configurations and different quantities 
of particles which, depending on the excitation energy 
of the system, can be evaporated before and after the 
breakup. 



2. Multifragmentation 

The multifragmentation mode extends over a large 
range of recoil velocities and indicates that the kinetic 
energy is shared among a large multiplicity of fragments 
with small mass asymmetry. In the analysis of the sys- 
tem 56 Fc(i J 4GcV) + P such a shape was proposed as a 
signature compatible with multifragmentation. In such 
system, which has a larger excitation energy per nucleon 
than 136 Xc+p, the multifragmentation mode was found 
to predominate already above carbon [Z(|. Along the 
same line, in agreement with the picture that the multi- 
fragmentation share increases with the excitation of the 
system, inclusive measurements at the incident energy of 
1 GeV per nucleon found only a small contribution from 
the Coulomb-hole component in comparison with the 
multifragmentation mode: it was the case of the systems 
56 Fe+ nat Ti [4O1 i36 Xe + nat Ti [§3, 136 Xe+ 208 Pb [&| , and 

124 Xe+ 208p£|gJ) 
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FIG. 10: (Color online) Evolution as a function of the mass 
number of the residue of the mean standard deviation of the 
invariant-cross-section distribution averaged on the atomic 
number. For lower masses than A — 30 the width associ- 
ated to the forward side of the multifragmentation compo- 
nent and the width of each peak of the Coulomb-hole mode 
are shown. The mean standard deviation of the full spectrum 
is shown for masses heavier than A = 30. The systematics 
of Morrissey [U and the prescription of Goldhaber [k| are 
shown for comparison. 



locity distribution confirms that the contribution of the 
Coulomb-hole component is still large. 



4- Nuclide production 



3. Widths of the velocity spectra 

In the system 136 Xe( lj 4GcV) + Pi the contribution of 
the multifragmentation mode prevails gradually on the 
Coulomb-hole mode for elements of larger mass number 
till, around silicon, the unfolding of the two components 
becomes too arduous. For silicon nuclides, the contribu- 
tion of the Coulomb-hole mode is still about half the pro- 
duction cross section of silicon and should still survive for 
a long range of elements. In this respect, the presence of 
the Coulomb-hole mode is still appreciable in the widths 
of the longitudinal- velocity spectra. In fig. [TUl the evolu- 
tion of the mean standard deviation (<xjp of the invariant- 
cross-section distribution is shown as a function of the 
mass number of the residue for A > 30. For A < 30 the 
mean standard deviation extracted from the forward side 
of the multifragmentation mode in the invariant-cross- 
section representation is plotted: this standard deviation 
is equal to the average standard deviation (crP ) of the 

II i ^ 

components of the multifragmentation mode. Around 
silicon, the shape of the multifragmentation mode loses 
its tail and the width extracted from the forward side 
coincides with the width of the full multifragmentation 
mode. The mismatch at A — 30 between the widths 
of the multifragmentation component and of the full ve- 



A survey on the production cross section, taken from 
ref. [HI, is presented in fig. QT] The contributions of 
the Coulomb-hole and multifragmentation modes, de- 
duced in this work, are indicated; a calculation is also 
included in the figure to estimate the contribution of a 
Weisskopf evaporation chain restricted to the emission of 
only neutrons, protons and alpha particles. Combining 
the production and the kinematic observables, we can es- 
timate that in the region of mass numbers corresponding 
to the minimum of cross section, while the information 
on the Coulomb-hole component gradually vanishes and 
the production cross section of intermediate-mass frag- 
ments drops steeply, the yield of the evaporation residues 
rises and the kinematics reflected in the velocity spectra 
is replaced by the characteristic Gaussian shape of the 
evaporation residues. Above half the mass number of the 
projectile we observe only the kinematics of the heaviest 
fragments. When approaching the mass of the projectile, 
as shown in fig. 1101 the velocity width of the invariant 
cross-section distributions evolves in qualitative agree- 
ment with the empirical systematics of Morrissey [53J, 
which is suited for describing the kinematics of heavy 
residues, and with the prescription of Goldhaber [52l|. 

The production cross section of evaporation residues 
(in this case we consider any evaporation pattern where 
intermediate mass fragments are not produced by fission 
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FIG. 11: (Color online) Projections of the residue-production 
cross sections along the mass number. For the light fragments 
the contribution from the two kinematic modes are indicated. 
A Weisskopf calculation for the expected contribution of the 
residues which evaporated neutrons, protons and alpha par- 
ticles is shown. 



in the first decay steps) varies monotonically as a func- 
tion of the mass number in agreement with the expec- 
tation that also the excitation energy of the hot frag- 
ments has a monotonic dependence (see fig. [T]). In this 
respect, the fact that the mass distribution of evaporation 
residues does not extend monotonically towards the light- 
est masses is an additional indication that intermediate- 
mass fragments are not fed by evaporation residues. We 
may remark that this conclusion, clearly suggested by 
the production cross sections, can not be inferred from 
the velocity spectra. There is in fact a gradual tran- 
sition between the shape of the velocity distributions 
of the heaviest intermediate-mass fragments, which are 
dominated by the multifragmentation mode and which 
are almost gaussian, and the gaussian-like spectra of the 
evaporation residues. For this reason, in lighter systems 
like 56 Fe+p [4(J, where intermediate- mass fragments and 
evaporation residues do not occupy distinguishable re- 
gions of the production cross-section distribution, wc 
could not identify the transition between the two pro- 
cesses from a similar experimental analysis. 

In ref. [45| a total reaction cross section of 1393 ± 
72 [mb] was obtained. We supposed that the Coulomb- 
hole mode corresponds essentially to binary splits. Al- 
though the production of the heavy fragments in the bi- 
nary splits could not be disentangled, it must be equal 



to the contribution of the corresponding light partners, 
which we measured for the greatest part, till A = 30, 
where the production cross section already drops to very 
small values. Under this assumption, the portion of total 
cross section related to the Coulomb-hole mode could be 
deduced: we found a value of 59.0 ± 3.7 [mb]. The Weis- 
skopf calculation of residues which evaporated protons, 
neutrons and alpha particles gives an indicative estima- 
tion of the corresponding fraction of total cross section 
of around 90%. The contribution of the multifragmenta- 
tion mode can not be deduced because neither the frag- 
ment multiplicity of the process, nor the total yield of 
this mode could be deduced from the measurement. The 
sum of the production cross sections related to the mul- 
tifragmentation mode which were measured up to silicon 
is of 8.6 ± 0.6 [mb]. By subtracting the measured con- 
tribution of the Coulomb-hole mode and the calculated 
fraction of evaporation residues to the measured total re- 
action cross section we estimate an indicative cross sec- 
tion of around 35/ (M) [mb] for the multifragmentation 
mode, where (M) is the corresponding mean fragment 
multiplicity, supposed to be larger than two. 



5. Fine structure in the nuclide cross sections 

A close analysis of the nuclide cross sections reveals a 
staggering of large magnitude for chains of nuclides hav- 
ing the same N — Z. This feature, which was already 
discussed in refs . l69l . [70l ] and analysed for the system 
136 Xe+p in ref. [45[ , is now presented in fig. [T^] for the 
two decay modes distinctly. The upper row presents the 
measured cross sections and the lower row quantifies the 




FIG. 12: (Color online) Upper row. Chains of nuclide cross 
sections with a given value of N — Z present an "even-odd" 
staggering. Lower row. Amplitude of the staggering analysed 
by the method of Tracy [Tl|] . It is positive when the staggering 
is "even-odd" and negative when the staggering reverses. 
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staggering amplitude according to the procedure intro- 
duced by Tracy [7l|. Despite the smoother evolution 
of the staggering strength for the Coulomb-hole mode, 
we observe the same general phenomenology for the two 
modes. In a simple Weisskopf picture this staggering 
was described as a manifestation of the very last steps of 
the decay path, independent of the entrance channel; the 
number of excited levels of the mother nucleus that could 
decay into the daughter nucleus determines the probabil- 
ity of a channel. In particular, even-mass nuclides show 
a smooth variation of the separation energy as a function 
of the element, once shifted by the pairing gap. There- 
fore, their formation is enhanced due to the pairing gap, 
which determines a larger number of excited levels of the 
odd-mass mother nuclei available for the decay into the 
ground states of even-even daughters rather than odd- 
odd daughters. This explains the even-odd staggering for 
N = Z and N = Z+2. On the other hand, the formation 
of odd-mass residues is not determined by the position 
of their ground states, which are ordered along the same 
mass parabola for one given mass number. Their forma- 
tion reflects the structure of the separation energy. Since 
neutron-rich nuclei have lower separation energy for neu- 
trons, their formation reflects the structure of the neu- 
tron separation energy and odd elements are favoured. 
This explains the reversal of the even-odd staggering for 
N = Z+l and N = Z+3. This simple picture is sufficient 
to conclude from the large amplitude of the staggering 
that the secondary decay affects largely and in the same 
manner the production of intermediate mass fragments 
from both decay modes. We disregard possible exotic de- 
cays which go beyond the Hauser-Fcshbach approach for 
compound-nucleus decay, like cluster emission (72l . p73| . 
and which might have an additional influence. 



B. Connection between the entrance channel and 
the exit channel 



The precise knowledge of the mean recoil velocity of 
the ending fragments of the decay process provides in- 
dications on the connection between the exit channel 
and the entrance channel. Such a study profits from the 
high-resolution measurement of the inclusive experimen- 
tal approach. However, this analysis can not be based 
on single-event observables, due to the large fluctuations 
induced by the Fermi momenta of abraded nuclcons and 
the recoil due to the decay. In order to smear out these 
effects, this analysis is based on employing, as an inclu- 
sive quantity, the mean value of the source velocity of one 
specific component of the velocity distribution in the pro- 
jectile frame. It is evident that a study of the entrance 
channel based on single-event observables, concerning for 
instance the information on the excitation energy, would 
go beyond the specificities of this inclusive experimental 
method. 



1. Mean recoil velocities 

Fig. [2] shows that a systematic displacement in recoil 
velocity between the two kinematic modes, fission and 
multifragmentation, could be resolved with high preci- 
sion since the two components are measured as parts of 
the same spectrum in the same set of experimental runs. 
The distinction of the two kinematic modes is a signature 
of the interplay of two distinct exit channels; their rela- 
tive displacement, which characterises all spectra of the 
intermediate-mass fragments shown in fig. [H attributes 
the two exit channels to different entrance channels. This 
displacement is in fact connected to the mean momentum 
transfer, which reflects the violence of the collision, and 
which evolves in different ways according to the decay 
mode. Such study is presented in fig. 1131 where three 
mass regions can be distinguished. 

1. Heavier masses than half the mass of the projectile 
(A > 68). They are necessarily the heaviest frag- 
ments produced in the reaction. They cannot be 
the lighter fragments of a binary decay, and they 
are not typical multifragmentation products, which 
should be more than two and should have compa- 
rable size. 

2. Intermediate-mass fragments with A < 30. Ex- 
plicit information on the source velocities of the two 
kinematic modes could be obtained experimentally. 
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FIG. 13: (Color online) Evolution as a function of the mass 
of the residue of the mean recoil velocity, deduced for the two 
kinematic modes separately for lower masses than A — 30. 
The mean value of the reconstructed full velocity distribu- 
tion is shown for the whole mass range. The systematics of 
Morrissey [kJ is shown for comparison. 



15 



3. Fragments in the range 30 < A < 68. This is in 
the minimum of the production as shown in fig. 111! 
only the average source velocity could be measured, 
without any kinematic distinction of the decay pro- 
cess. 

The region of heavier mass numbers than half the mass 
of the projectile is dominated by the presence of evap- 
oration residues. Here, the mean recoil v?, evolves al- 
most linearly with the mass loss from the projectile ve- 
locity down to around —0.2 cm/ns with respect to the 
projectile frame. The intcrmcdiatc-mass production at- 
tributed to the Coulomb-hole mode should originate from 
the decay of heavy residues which evaporate not only 
nucleons but also heavier fragments. In the process, 
the accompanying heavy residues should have suffered 
a large loss of mass. Therefore, it is reasonable that the 
largest negative mean source velocities, of around —0.25 
cm/ns, reflected by intermediate-mass fragments issued 
of binary splits, are close to the source velocities which 
characterise residues of around half the projectile mass 
(A ~ 68). With less statistical significance, the data 
on intermediate-mass fragments from the coulomb-hole 
mode suggest also the tendency of reducing the mean re- 
coil with increasing mass: the mean recoil varies from 
around —0.2 cm/ns for A ~ 5 to around —0.25 cm/ns 
for A ~ 30, in the projectile frame. This is reasonable, 
because the emission of heavier intermediate-mass frag- 
ments from binary decay is expected for higher excitation 
energies; at the same time, higher excitation energies cor- 
respond to smaller impact parameters, reflected by less 
heavy partners which also tend to have larger negative 
mean velocities. The mean recoil of the intermediate- 
mass fragments issued of multifragmcntation could be 
measured in a limited mass range (A < 30). For the 
heaviest fragments, it approaches almost the velocity of 
the projectile around A = 30. For the lighter fragments, 
we find a strong tendency to larger negative values in the 
projectile frame. Both the rather high velocities of mul- 
tifragmcntation products around A = 30 and the strong 
trend to very large negative values towards A ~ 10 have 
no reasonable counterpart in the heavy residues. This is 
not surprising, since, on the one hand, the heavy frag- 
ments cannot be produced in a typical multifragmcnta- 
tion process; on the other hand, also the convex pat- 
tern in the distribution of invariant cross section of the 
multifragmcntation products contradicts the simultane- 
ous emission of a heavy partner. For larger masses than 
around A = 15, the source velocities of the intermediate- 
mass fragments from multifragmcntation variate in the 
same range of the source velocities of both the heavier 
and the lighter partners of binary splits; still, the two 
decay modes can be distinguished by the shape of the 
velocity distribution. 

Altogether, the experimental results on binary-decay 
events suggest that the source velocities can be well ex- 
plained by the empirical systematics of Morrissey [53| : 
both in the heavy (A > 68) and in the light (Coulomb- 
hole mode measured up to A = 30) residues, there is a 



tendency towards more violent collisions from the heav- 
iest, respectively lightest, masses towards more symme- 
try. In particular, in the case of binary decays, three 
quantities are correlated: the evolution of the system to- 
wards more symmetric masses is a signal of increasing 
friction for less peripheral impact parameters, and is re- 
flected in the mean-source velocity. In the case of light 
intermediate-mass fragments issued of multifragmcnta- 
tion, the correlation between two of these quantities, the 
size and the impact parameter, is also well known from 
previous experiments [74j , where the fragment size was 
observed to decreases with increasing energy deposition, 
for more central collisions. 



2. The mean recoil velocities as a possible indicator of the 
impact parameter 

The new additional information of this experiment, 
that the size of the intermediate-mass fragments corre- 
lates with the mean recoil in the whole mass range and for 
both the two kinematic modes, makes plausible to estab- 
lish a further correlation between the mean source veloc- 
ity and the impact parameter also for the multifragmen- 
tation products. We think this is an important generali- 
sation of the interpretation of the kinematic information, 
which combines results from previous experiments with 
the new findings of this work. In this respect, this experi- 
ment extends the well known Morrissey kinematics which 
was established for masses very close to the heavy reac- 
tion partner, all the way down to the intermediate-mass 
fragment region. Using the well established correlation of 
the fragment size in the multifragmcntation regime with 
the impact parameter, we show that the fragment size 
also correlates with the source velocity. This leads to the 
conclusion that the source velocity is another variable 
suitable to characterise the impact parameter. 

However, differently from the case of the heavy 
residues, the absolute mean source velocity tagged by 
the intcrmcdiatc-mass fragments produced in multifrag- 
mentation may manifest a reversed dependence, reducing 
to very small values in the projectile frame. The exper- 
imental disentangling of the two kinetic modes could be 
extended up to rather heavy fragments: it allowed to 
make evident the presence of a reversal point of such de- 
pendence in the region of the heaviest intermediate mass 
fragments, around A = 25, determined by the multifrag- 
mcntation products. Such finding results from the high- 
resolution of the inclusive method. For heavier masses, 
the measurement could no more distinguish the two ki- 
netic modes. However, the reversal is also evident, even 
if smoothed, by analysing the overall behaviour of the 
sum of the two kinetic modes, which could be measured 
for the full range of masses; such observable connects 
the evolution of the intermediate-mass fragments with 
the behaviour of the heavy residues. Due to this ef- 
fect, the impact parameter is not uniquely defined by 
the mean source velocity in the projectile frame for any 
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fragment size. This ambiguity, can be controlled by ad- 
ditionally accounting for the information on the size of 
the largest fragment: the products which appear with 
the same mean velocity in the projectile frame are very 
much separated, if we consider the mass of the heaviest 
fragment formed in the reaction as an additional crite- 
rion. In this respect, a correlation of the mean velocity 
in the projectile frame with the impact parameter could 
be established also in the intermediate-mass region. 



3. Possible traces of phase transitions? 

Supposing that in proton-induced reactions at high in- 
cident energies dynamical effects arc negligible, we find 
that the dependence of the recoil observable as a function 
of the mass number is very different for the Coulomb-hole 
and multifragmcntation modes. As depicted in fig. [HI for 
the first mode it varies in a reduced range, while for the 
second it presents a strong dependence. 

Keeping in mind that the recoil observable can be in- 
terpreted as an indicator of the impact parameter, this 
behaviour may be interpreted tentatively as follows. The 
Coulomb-hole mode is associated to a small range (with 
respect to the uncertainties) of recoil velocities in the pro- 
jectile frame: this may indicate that, within the picture 
of a Liquid-Fission phase transition, this mode is also as- 
sociated to small energy fluctuations when the system is 
situated in the coexistence region, and it corresponds to a 
phase transition with vanishing latent heat, as expected 
for the Liquid-Fission phase transition [HI, H3, HH, IH| - 
The situation changes completely for the multifragmen- 
tation mode, which is associated to a very large range of 
recoil velocities in the projectile frame; this may corre- 
spond to a strong dependence of the fragment configu- 
ration on the excitation energy, which corresponds to a 
finite latent heat. If we consider the heavy residues and 
the intermediate-mass fragments issued of multifragmen- 
tation, and we exclude the Coulomb-hole mode, the plot 
of fig. [T3] gives a good information on the evolution of 
the recoil as a function of the heaviest fragment. With 
increasing absolute values of the average recoil in the pro- 
jectile frame, the system firstly explores a pattern charac- 
terised by only heavy evaporation residues, close to the 
mass of the projectile; it then explores a configuration 
where both the intermediate-mass-fragment production 
and the residues coexist; finally, it evolves towards a situ- 
ation where the very large average absolute recoil velocity 
in the projectile frame corresponds only to the multifrag- 
mcntation mode and no more to the residue production. 
These inclusive observables are not incompatible with 
signals of bimodality, which were widely investi gate d in 
recent exclusive experiments [IjJ [2(| [23, HE US US HH, 
and which are proposed as a manifestation of the liquid- 
gas phase transition in nuclear systems. Our interpre- 
tation, based on average observables, follows closely this 
line, but it should be corroborated and quantified by ex- 
clusive experiments. 



VI. CONCLUSIONS AND PROSPECTS 

In this report we developed a specific procedure for 
the deconvolution of the distributions of longitudinal ve- 
locities of intermediate-mass fragments, which were mea- 
sured inclusively and identified in atomic and mass num- 
ber. Technically, the experimental approach we adopted 
in measuring the characteristics of the intcrmcdiatc-mass 
fragments has several unique advantages: it avoids the 
contamination from target-like reaction partners with 
low kinetic energy which is a difficulty in direct kine- 
matics experiments, and the overloading by the primary 
beam at small forward angles which challenge inverse- 
kinematics experiments in absence of a spectrometer; it 
also profits from the absence of the low-energy detection 
threshold which characterises silicon detectors, and prof- 
its from the absence of magnetic aberrations of the achro- 
matic setting of the small-acceptance spectrometer. As a 
main disadvantage, the reduced acceptance precludes to 
measure particle correlations. 

As a mayor general result, we could disentangle and 
analyse quantitatively, in terms of invariant cross sections 
and in terms of production cross sections two kinematic 
modes for the exit channel in the spallation of 136 Xe at 
1 GeV per nucleon, for several isotopes of light elements 
ranging from lithium to silicon. 

One mode becomes manifest in a Coulomb shell in ve- 
locity space and a convex shape in the distribution of 
invariant cross section, related to events with low frag- 
ment multiplicity and large asymmetry in mass number. 
This mode is attributed to asymmetric binary decays, 
and the fission velocity deduced from the measurement 
is consistent with the binary split of a system close to 
the mass of the projectile; this is however less true for 
elements lighter than carbon, which exhibit very high 
velocities, incompatible with standard calculations. The 
results are in line with previous similar studies of the 
formation of intermediate-mass fragments in the lighter 



system 56 Fe+p [40( and in the heavier system 38 U+p [61 1 



at the same incident energy of 1 GcV per nucleon. 

Another mode manifests with a convex shape in the 
distribution of invariant cross section, related to events 
with large fragment multiplicity and small asymmetry 
in mass number. This mode is attributed to multifrag- 
mentation. As a relevant result of the present analysis 
we indicate that, despite the small fraction of the total 
cross section which feeds multifragmcntation, the contri- 
bution of this process to single nuclide cross sections is 
found to be large for intermediate-mass fragments and 
even dominant for some specific nuclides. Such a large 
contribution and, more generally, the large production 
of light elements, may have consequences for technologi- 
cal applications of spallation, such as transmutation and 
safety. 

The main result of this work is the identification of 
these two decay processes, and the extraction of kine- 
matic and production observables related to the exit 
channel and associated to the two processes separately 
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We measured the characteristics of the ending fragments 
of the decay process, which give valuable information on 
the heaviest reaction product; the latter is even a funda- 
mental probe of the thermodynamic properties of finite 
nuclear systems, we also concentrated on the sequence 
of evaporation decays, reflected in the staggering of nu- 
clide cross sections; such staggering has a comparable 
amplitude for the two decay modes, and we attributed 
it to a large fraction of secondary decay; unfortunately, 
such process has the effect of smearing out the observ- 
ables which could connect the entrance channel to the 
exit channel. 

However, signals of this connection are found back in 
the dependence of the mean recoil as a function of the 
mass number. By combining results from previous ex- 
periments with the new information of this work, we es- 
tablished the inclusive mean source velocity valu- 
able indicator of the impact parameter. By relying on 
such indicator, we firstly confirmed that the Coulomb- 
hole and the multifragmentation modes, reflect different 
processes. Then, on the more speculative side, we placed 
our study in the context of the recent experimental and 
theoretical investigations on the phase transitions in nu- 
clear systems. We advanced the idea that the two modes 
could also be related to the phase-transition phenomenol- 
ogy in the finite nuclear system. The Coulomb-hole 
mode exhibits compatible features with a Liquid-Fission 
phase transition, characterised by vanishing latent heat. 
The characteristics which describe the multifragmenta- 
tion mode are in qualitative accordance with recent ob- 
servations of a bimodal distributions of the largest frag- 
ment charge: for the less violent collisions the multifrag- 
mentation mode does not appear and only the heaviest 
evaporation residues are produced. The multifragmen- 
tation mode then gradually prevails for more and more 
violent collisions till it dominates the decay pattern for 
the most violent collisions, leading to the production of 
the smallest residues. This behaviour was found in the 
study of projectile-like fragments formed in (semi-)pe- 
ripheral collisions of heavy ions at Fermi energies and it 
was based on exclusive event selections; it was interpreted 
as a robust signal of the Liquid-Gas phase transition in 
the finite nuclear system [1^, [H, H3, [H, IM IM HI • In 
this case, the dense phase was identified with residue-like 
fragments, after the suppression of fission events. 

Relativistic nuclear reactions induced by hadrons were 
already employed for investigating the statistical nature 
of the decay process [ll|. The fraction of incident en- 
ergy which is transferred to the prefragmcnt appears al- 
most completely as single-particle excitation, while col- 
lective effects like compression are relatively small; for 
this reason, we suggest that spallation reactions may 
be also suited for distinguishing the signals of Liquid- 
Fission and the signals of the Liquid-Gas phase transi- 
tions in finite nuclear systems. In the same context, we 
advance the view that a method of inclusive event selec- 
tion like the one presented in this paper and, more gen- 
erally, event selections based on kinematic observablcs, 



could be relevant in a further generation of experiments 
also at Fermi energies, where the measurement of iso- 
topic and kinematic observables in correlation with light- 
fragments is a mayor request and a foreseen development 
of the previous-generation exclusive experiments. This 
is already the direction followed by very recent experi- 
ments, where a spectrometer is coupled to a Att detec- 
tor array [75j . In this respect, high-precision momen- 
tum measurements and the nuclide identification of the 
residues, from inclusive approaches, can offer a further 
tool of event selection if coupled to correlation measure- 
ments, from exclusive approaches. 



APPENDIX A: VELOCITY RECONSTRUCTION 
WITH MULTIPLE SOURCES 

In this section we give a general description of the de- 
convolution of the measured longitudinal velocity distri- 
butions shown in fig. [21 which leads to the invariant cross 
sections shown in fig. 0J This procedure is based on the 
three conditions listed in section IIIII For generality, in 
this section we suppose that one velocity spectrum re- 
sults from a continuous distribution of sources along the 
beam axis. 

The deconvolution imposes to change among three ref- 
erence frames, the laboratory frame, the projectile frame, 
and the reference frame of one given isotropic source, la- 
belled L,P,u, respectively. We define the velocity vectors 
v p , u p , and v u = (v — u), to indicate the velocity of 
fragments in the projectile frame, the velocity of a given 
source u in the projectile frame, and the velocity of frag- 
ments in the frame of a given source u, respectively 

We intend to pass by deconvolution from the measured 
longitudinal velocity distribution d/(t;jp/di>jj', which de- 
pends on the acceptance, to the (non-relativistic) invari- 
ant cross section ai(v p ), which is not depending on the 
acceptance. Since we must proceed through unfolding 
over the distribution of sources, we define the correspond- 
ing distributions associated to a given source u as in the 
following scheme: 

full space source u 
acceptance dl{v\)/dv\ d[8 u I{vfj\/dv§ 
no acceptance cr\(v v ) 8 u a\{v u ) 

d[<5 qj ,-T('u")] / , where 5 U stands for d/dit p , is the portion 
of the measured longitudinal velocity distribution which 
we attribute to a given source u, and 5 u a\{v u ) is the 
portion of the invariant cross section which we attribute 
to the source u. 

The measured production yield d/(vy > )/d?;jj > of one indi- 
vidual nuclide is the integration over the emitting sources 
of the quantity 5 u I(v p ), which is the portion of the 
spectrum associated to the sources lying in the interval 
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[ U p,uP + dv?]: 



^1(1^) d«P. 



(Al) 



The angular acceptance of the magnetic spectrometer, 
with an almost circular aperture of ((^(ip)) » 15 mr 
in the laboratory frame (the azimuthal angle ip around 
the beam axis is introduced to take into account the de- 
viation of the acceptance shape from a circle) imposes 
that a limited selection of the emission-velocity space 
could be explored and only the selected part of 5 U I ap- 
pears in the measured yields. The Lorentz transforma- 
tion attributes to this selected portion a conical bound- 
ary along the beam axis in the frame of one emitting 
source, defined by the velocity vector vj_, of magnitude 
Vj_ = 7^(w L + v^)tga L (<p), and orientation along the az- 
imuthal angle p> in the plane orthogonal to the beam axis; 



is the longitudinal projection of the boost v u in the 
frame of the source u, imparted to the emitted fragment. 
By moving to the frame of the source u, the contribution 
d[(5 u /(wp]/dt'^ of the source u to the measured yields is 
described by the integration of the transverse velocity Vj_ 
over the portion of the plane selected by the ensemble of 
vectors so that 



dvl? 



2tt I «5. 1 



x dv u x 



(A2) 



where a is the production cross section for the nuclide 
associated to the velocity spectrum. Since in our case 
the velocities in the source frame are not relativistic, the 
argument of the integral in cq. (|A2[) is equivalent to the 
invariant cross section as a function of the velocity in the 
source frame: 



d[S u a(v u )} 
dv u 



w? EU d[5 u a(v")] 



dp u 



-J„ai(« u ), (A3) 



where E u and p u indicate the total energy and the mo- 
mentum in the source frame, respectively. Introducing 
the invariant cross section in the integral we write: 



d[SuI(vtt)} 



2tt \vl\ 



dw™ 



6 u <n(v u )vldvldtp, (A4) 



This expression does not impose any constraint for 
the distribution of emitted fragments. Under the 
assumption of isotropy S u ai(v u ) = 5 u <Ji(v u ) = 
[l/Air(v u ) 2 ]d[S u a(v u )]/dv u and by substituting the inte- 
gration variable u" with v u = \J( V ™) 2 + ( v ± ) 2 we can 
write (indicating the vector of magnitude v^, oriented 
along the beam direction) : 



where the integral boundaries are \vl 



i p | and 



- MP) 2 + {vl) 2 . 

The distribution of invariant cross section ai(v u ), pre- 
sented in fig. [4] was deduced by the deconvolution of the 
integral eq. (|A5|) . The deconvolution can be solved only 
by restricting to a finite number of degrees of freedom: 
this implies that the number of sources is limited (the 
integral over u p is converted into a discrete sum over the 
sources). With this general prescription we constructed 
an axial-symmetric distribution of invariant cross section 
in the beam frame cri(t; p ) as the superposition of isotropic 
emission patterns, each one defined by one individual 
source: 



da 



i 2 dvP 



= a:K) = J 5 u a Iu (\vP-v?\)dv?. 



(A6) 



This spectrum would be obtained directly from the ex- 
perimental measurement only if two conditions were sat- 
isfied. First, the angular acceptance should not impose 
any cut. Second, also the transversal components of the 
velocity should be measured. 



APPENDIX B: PARAMETERS OF THE 
INVARIANT CROSS SECTIONS 

The analysis described in this work results from the 
deconvolution procedure introduced in section IIIII and 
applied to the analysis of section lTVl so as to produce the 
results sketched in fig. [3l As an outcome, the procedure 
gave the parameters listed in table IIIII 



dl(v^ 
dv* 
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5 u ai(v u )v u dv u dtpdu p , (A5) 
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TABLE III: Parameters of the invariant cross sections, which determine the plots in fig.|4j fig. HOI and fig. 1131 The parameters 
are the mean recoil of the Coulomb-hole mode , the ridge of the ring associated to the Coulomb-hole mode in the invariant 
cross section distribution f poa k, the standard deviation of the Coulomb- hole mode in the invariant cross section distribution 
Cm^i the mean recoil of the multifragmentation mode (u^), and the average standard deviation of the components of the 
multifragmentation mode in the invariant cross section distribution (erf ). 
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